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Abstract. We present new grids of pre-main sequence 
(PMS) tracks for stars in the mass range 0.1 to 7.0 M©. 
The computations were performed for four different metal- 
licities (.2=0.01, 0.02, 0.03 and 0.04). A fifth table has 
been computed for the solar composition (2=0.02), in- 
cluding a moderate overshooting. We describe the update 
in the physics of the Grenoble stellar evolution code which 
concerns mostly changes in the equation of state (EOS) 
adopting the formalism proposed by Pols et al. (1995) and 
in the treatment of the boundary condition. Comparisons 
of our models with other grids demonstrate the validity of 
this EOS in the domain of very low-mass stars. Finally, we 
present a new server dedicated to PMS stellar evolution 
which allows the determination of stellar parameters from 
observational data, the calculation of isochrones, the re- 
trieval of evolutionary files and the possibility to generate 
graphic outputs. 



"i> 1. Introduction 



Interests in modeling of very low-mass stars (VLMS) has 
grown steadily over the past decade due to big efforts in 
both theory and observations. Brown dwarfs and very low- 
mass stars are now frequently discovered and the abun- 
dant literature in this domain attests of the vitality of 
this field. From these observations, we now realize that 
very low-mass stars represent a large fraction of the stel- 
lar population. However, the precise determination of the 
initial mass function, and more generally of stellar param- 
eters and evolutionary states is still limited by the accu- 
racy of stellar models. It is therefore desirable that differ- 
ent tracks be available in order to understand the theoret- 
ical discrepancies and also to estimate the uncertainties 
associated with the derivation of stellar parameters. 

On the theoretical side, the computation of the struc- 
ture of very low-mass stars is very challenging since it in- 
volves the micro-physics associated with dense, cold and 
partially degenerate matter. In the interior of VLMS, col- 
lective effects due to large densities become important and 
an accurate treatment of pressure ionization and Coulomb 
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interactions is required. Also, in the cool atmosphere of 
VLMS, several molecules form and produce strong absorp- 
tion bands that considerably modify the emergent spec- 
trum. As a consequence, the stellar surface cannot be con- 
sidered as a black body and stellar atmosphere models 
must also be used. 

Since our last grids (Siess et al. 1997, hereafter SFD97), 
substantial improvements have been made in the Grenoble 
stellar evolution code that led us to the production of new 
grids of PMS tracks. The main modification to the code 
concerns the equation of state (EOS) which treats more 
accurately the behavior of cold and degenerate matter. In 
particular, we now take into account the effects associated 
with the pressure ionization and treat more carefully de- 
generacy conditions in stellar interior. We also updated 
the opacity tables using the latest release of the OPAL 
group and modified our surface boundary conditions using 
analytic fits to stellar atmosphere models. The updates to 
the code are describe in Sect. ^ followed by a brief descrip- 
tion of the grids. Then in Sect. ^ we compare our models 
to grids recently published by other groups. Finally, we 
describe a new server dedicated to PMS evolution. 

2. The stellar evolution code 

A general description of the stellar evolution code has al- 
ready been presented in several papers (i.e. Forestini 1994, 
SFD97) and only subsequent modifications are discussed. 
They include update in the opacity tables, in the treat- 
ments of the atmosphere and the incorporation of a new 
formalism for the equation of state. 

2.1. Opacity tables 

At high temperature (i.e. above 8000K), we use the OPAL 
radiative opacities by Iglesias and Rogers (1996) which 
now include additional metal elements. The tables are in- 
terpolated in temperature, density and H, He, C, O and 
metal mass fraction using the subroutine provided by the 
OPAL group. However, for the computation of very low- 
mass stars (M ^0.2 M©), we had to extrapolate the opac- 
ity tables to logi? = \og{p/T%) = 4.0 but this is of little 
effect since the internal region concerned with this extrap- 
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olation (T 10 5 K) has an almost adiabatic behavior. 
Consequently, the temperature gradient is governed by the 
thermodynamics coming from the equation of state and 
is nearly independent of the opacity coefficient. Conduc- 
tive opacities are computed from a modified version of the 
Iben (1975) fits to the Hubbard and Lampe (1969) tables 
for non-relativistic electrons, from Itoh et al. (1983) and 
Mitake et al. (1984) for relativistic electrons and from for- 
mulae of Itoh et al. (1984) as well as Raikh and Yakovlov 
(1982) for solid plasmas. Below 8000K, we use the atomic 
and molecular opacities provided by Alexander and Fer- 
gusson (1994) as in SFD97. 

2. 2. Treatment of the atmosphere 

The stellar structure is integrated from the center to a 
very low optical depth (t = 0.005) in the atmosphere. 
In the regions where r < 10, we constrain the atmo- 
spheric temperature profile (as well as the radiative pres- 
sure and gradient) to correspond to those coming from 
realistic atmosphere models determined by the integra- 
tion of the radiative transfer equation. More specifically, 
we have performed one analytic fit of T(t), func- 
tion of effective temperature T e g, surface gravity g e g and 
metallicity Z which agrees with the atmosphere models 
to within 20% or less. The parameters of this fit are con- 
strained by various atmosphere models, namely (1) Plcz 
(1992) for 2500 < T cS < 4000K, (2) Eriksson (1994, pri- 
vate communication, using a physics similar to Bell et al. 
1976) up to 5500K and (3) Kurucz's models (1991) com- 
puted with the ATLAS 12 code above 5500K. Unfortu- 
nately, the lack of atmosphere models in some regions of 
the Z — g e g plane did not allow us to precisely determine 
the behavior of this fit at high gravity and/or non solar 
metallicities. This can potentially affect the modeling of 
very low-mass stars but considering the weak dependence 
of our fit on Z and g e s compared to the one on T e g, we de- 
cided to use this unique fit in all our computations instead 
of using a grey atmosphere approximation. This fit also 
provides a faster convergence in the surface layers due to 
the smooth profile of T(t) and of its derivatives. Finally, 
let us mention that we define the effective temperature 
as T e ff = (L/AttuR 2 ) 1 / 4 where L and R are estimated at 
t = 2/3. 

2. 3. Equation of state 

In order to improve the description of the micro-physics 
relevant to VLMS, we have revisited the computation of 
the equation of state (EOS), adopting the scheme devel- 
oped by Pols et al. (1995, hereafter PTEH). This new 
formalism is characterized by three main points : 

— we adopt the chemical point of view (i.e. Saumon et al. 
1995) where bound configurations (ions, molecules) are 
populated according to matter-photons and matter- 
matter interactions, 



— we use the assumption of local thermodynamic equi- 
librium which allows the use of the Saha equations to 
determine the ion abundances as a function of the lo- 
cal values of temperature, T and number of electrons 
per unit mass, N e , 

— the plasma is described by three components, photons, 
ions and electrons, coupled only through the photo- 
ionization and photo-dissociation processes. Each com- 
ponent of the plasma is represented by an independent 
term in the Helmoltz free energy (Fontaine et al. 1977). 
Non-ideal corrections due to Coulomb shielding and 
pressure ionization are treated separately in an addi- 
tional term which expression has been fitted to the 
MHD EOS (Mihalas et al. 1988). This fit allows us to 
treat the regime of strong coupling while keeping the 
formalism of a separate additive term in the free energy 
to describe non-ideal effects (as explained in Saumon 
et al. 2000). 

This new equation of state allows us to better take into 
account the effects of pressure ionization, partially degen- 
erate matter and Coulomb interactions and provides very 
smooth profiles of the thermodynamic quantities since an- 
alytical expressions of all their derivatives can be obtained. 
Note that, as we will see in the next section, the computa- 
tion of this EOS naturally introduces a new independent 
variable related to degeneracy. 

2.3.1. A new variable 

The general philosophy of this new EOS is to replace the 
density p (or pressure P) by a new independent variable 
which describes the electron degeneracy. To do so, we use 
the electron degeneracy parameter n e — /3£ e where [3 = 
1/kT and £ e is the electron chemical potential. In this 
context, the density is given by 

9 N e ( Ve ,T,y x y { > 

where n e , N e and yx are the number of electrons per unit 
volume, per unit mass and the molar fraction of a chemical 
element X, respectively. 

N e depends on n e and T through the generalized Saha 
equations, and through the mass and charge conservation 
equations 
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example), x l+ its ionization potential, Z xi + its partition 
function and Ax its atomic molar mass. Ari e contains the 
contributions owing to the non-ideal effects of Coulomb in- 
teractions and pressure ionization which will be described 
in the following section. The use of rj e instead of p as input 
variable for the Saha equations allows us to analytically 
derive the ionization states for all the desired elements 
anywhere in the star. 

The number of electrons per unit volume n e is then 
obtained in terms of r\ e and T through the Fermi-Dirac 
integral 

Stt [°° y 2 dy 

n - = aT / x^Z — ( ? ) 
A c Jo Ve + 1 

where y = = x 2 + 2x, A c is the electron Compton 
wavelength and T*—kT/m e c 2 . This is a general definition 
which allows for rclativistic electrons. In order to speed 
up the computation of n e , we use polynomial fits to this 
integral. These fits have a simple form when expressed in 
terms of two variables, called / and g by Eggleton et al. 
(1973), arising from the asymptotic expansion of Eq. ^. 
Following these authors, 

Ve = ln(/) + 2 x (v/T+7 - ln(l + vT + 7))- (8) 

As a consequence, the variable ln(/) (hereafter In/), very 
well suited to describe the degeneracy of the plasma, is 
our new independent variable replacing p in the stellar 
structure equations. The change of variable p(\nf, T, y{) 
is defined by Eqs. [fl, 0, ^ and ||. In weakly degenerate con- 
ditions, In/ can be identified with a good approximation 
to the degeneracy parameter r\ e while at high degeneracy 
/ scales as r\ e . 

2.3.2. Physical ingredients of the EOS 

The plasma is always considered as partially ionized for H, 
He, C, N, O and Ne. The other species are considered to be 
either neutral or totally ionized. Our tests indicate that 
for M < 3 M , the structure and evolution of the stars 
are almost indistinguishable in both prescriptions (as far 
as PMS is concerned). In the more massive stars (M ^ 
3 Mo) we report that models assuming total ionization 
for the heavy elements have initially a slightly shallower 
convective region but when the star reaches the ZAMS, 
these differences vanish. 

Non-ideal effects include Coulomb shielding for all spe- 
cies and pressure ionization for H, He, C, N and O, ex- 
cept for H . Both effects have been incorporated following 
PTEH by means of analytic fits to the non-ideal terms in 
the Hclmoltz free energy. The fit to the pressure ioniza- 
tion has been derived from the MHD EOS for a mixture 
AT=0.70 and Z=0.02. Nevertheless, they have been used 
in all our computations, regardless of the chemical com- 
position. According to Pols (private communication) , this 
is a reasonable approximation since these fits give the ex- 
pected behavior of pressure ionization as a function of T 
and of the degeneracy parameters (rj e or In/) . 
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Temperature from surface to center in Kelvin 

Fig. 1: Profiles of H, H+, H~ and H 2 inside a 0.1 M Q on 
the ZAMS. The star has an age of 8.36 x 10 9 yr, a solar 
composition and is located in the HRD at T cr r = 2776 K 
and L = 8.2 x 1CT 4 L 

The computation of H 2 and H~ abundance is taken 
into account. For H 2 , the partition function of Irwin (1987) 
is used in Eq. [?] below 16000K and had to be extrapolated 
at higher temperatures to ensure the convergence of the 
model. This is certainly of little importance since in this 
region H 2 is insignificant. It must be noted that for very 
low-mass stars (M ^ 0.4 M Q ), H 2 becomes very abundant 
in the surface layers and must be included in the com- 
putations to ensure consistency between the atmosphere 
models (which include this molecule) and the EOS. For 
example, in a 0.1 M Q 50% of the hydrogen is trapped into 
H 2 at a temperature of 11800K (Fig. [l]). Further inside 
the star, this molecule cannot survive. 
Due to its negative charge, the abundance of H _ derived 
from Eq. ||is proportional to exp(?y e ) instead of exp(— r] e ). 
This different dependence of H~ abundance on the de- 
generacy parameter is not reproduced by the fit to the 
pressure ionization, thus preventing for the moment the 
treatment of this effect on this ion. This different behavior 
also explains why the peak of H _ is located at the maxi- 
mum of the degeneracy parameter, around T = 5 x 10 4 K 
where T) e = 3.02. We also verified that because of its small 
abundance in the encountered conditions, H has a negli- 
gible effect on the evolution tracks, even if its influence on 
the opacity coefficient can be large. Finally, in the central 
region of the star, H only subsists in its ionization state 
H+. 

Compared to Pols et al. (1995), our EOS presents slight 
modifications which mostly concerns the treatment of ion- 
ization. More specifically, now we compute the ionization 
states of C, N, O and Ne everywhere in the star, the ele- 
ments for which ionization is not followed can be consid- 
ered either completely ionized or neutral, the most recent 
partition function of Irwin (1987) is used for H 2 instead 
of Vardya's (1960) and the presence of H~ has been ac- 
counted for in the Saha equations. Apart from these mod- 
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ifications, the basic scheme and physical inputs (mainly 
the treatments of the non-ideal corrections and degener- 
acy) are the same. We checked that by adopting the same 
ionization treatment as PTEH (i.e. without the computa- 
tion of H - and assuming complete ionization for all the 
elements heavier than He), all the physical quantities de- 
rived form our EOS are comparable to the one computed 
by PTEH to within less than 1% . Consequently, we refer 
the reader to these authors for a detailed comparison of 
this EOS with other ones. In particular, they found that 
due to the use of a simplified treatment of pressure ioniza- 
tion, this EOS can differ by up to 10% with the MHD and 
OPAL EOS, and even a bit more in the region of the p — T 
plane delimited by the 3.0 < logT < 4.5 and < \ogp < 2. 

3. Pre-Main sequence computations 

3.1. Initial models 

Our initial models are polytropic stars that have already 
converged once with the code. The central temperature of 
all our initial models is below 10 6 K so deuterium burning 
has not taken place yet. The stars are completely con- 
vective except for M > 6M Q where a radiative core is 
already present. For solar metallicity (Z — 0.02), we use 
the Grevesse and Noels (1993) metal distribution. For dif- 
ferent Z, we scale the abundances of the heavy elements 
in such a way that their relative abundances are the same 
as in the solar mixture. 

3.2. The grids 

Our extended grids of models includes 29 mass tracks 
spanning the mass range 0.1 to 7.0 M Q . Five grids were 
computed for four different metallicities encompassing 
most of the observed galactic clusters (Z = 0.01, 0.02, 
0.03 and 0.04) and, for the solar composition (Z = 0.02), 
we also computed a grid with overshooting, characterized 
by d = 0.2H p , where d represents the distance (in unit of 
the pressure scale height H p measured at the boundary 
of the convective region) over which the convective re- 
gion is artificially extended. An illustration of these grids 
is presented in Fig. ||. Our computations are standard in 
the sense that they include neither rotation nor accretion. 
The Schwarzschild criterion for convection is used to de- 
limit the convective boundaries and we assume instanta- 
neous mixing inside each convective zone at each iteration 
during the convergence process. 

During the PMS phase, the completely convective star 
contracts along its Hayashi track until it develops a radia- 
tive core and finally, at central temperatures of the order 
of 10 7 K, H burning ignites in its center. The destruction of 
light elements such as 2 H, 7 Li and 9 Be also occurs during 
this evolutionary phase. Deuterium is the first element to 
be destroyed at temperatures of the order of 10 6 K. The 
nuclear energy release through the reaction 2 H(p,7) 3 He 
temporarily slows down the contraction of the star. Then, 
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Fig. 2: Evolutionary tracks from 7.0 to 0.1 M for a solar 
metallicity (Z=0.02 and y=0.28). Isochrones correspond- 
ing to 10 6 , 10 7 and 10 s (dashed lines) arc also represented. 
This figure has been generated using our internet server 



7 Li is burnt at ~ 3 x 10 6 K shortly followed by 9 Be. For 
a solar mixture, our models indicate that stars with mass 
< 0.4 M Q remain completely convective during all their 
evolution. We also report that in the absence of mixing 
mechanisms, stars with M > 1.1 never burn more 
than 30% of their initial 7 Li. 

The account for a moderate overshooting character- 
ized by d = 0.20H p significantly increases the duration of 
the main sequence (MS) for stars possessing a convective 
core and provides additional Li depletion during the PMS 
phase. More quantitatively, the MS lifetime of a 1.2 M Q 
star is increased by 25%, this percentage then decreases 
to level off around a 15% increase for M 3M Q . Surface 
depletion of Li occurs only in stars with 0.4 < M ^ 1.3 
for which the temperature at the base of the convective 
envelope can reach 3 x 10 6 K. Our models including over- 
shooting indicates that Li is much more efficiently de- 
pleted around 0.8 M Q , where its abundance is ~ 3 x 10 5 
smaller than in a standard evolution. 

3.3. Our fitted solar model 

With the above physics, we have fitted the solar radius, 
luminosity and effective temperature to better than 0.1% 
with the MLT parameter a — 1.605 and an initial com- 
position Y = 0.279 and Z/X = 0.0249 that is compatible 
with observations. The values used for this fit are quite 
similar to those obtained by other modern stellar evolu- 
tion codes (see e.g. Brun et al. 1998). Our sun, computed 
in the standard way, i.e. without the inclusion of any dif- 
fusion processes, has the following internal features 

- at the center : T c = 1.552 x 10 7 K, p c = 145.7 gem" 3 , 
Y c = 0.6187 and the degeneracy parameter r\ c = 
-1.527; 
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Table 1: PMS stellar models used in the comparisons 



X 


Z 


group 


code 


0.725 


0.019 


Baraffe et al. (1998) 


BCAH 


0.680 


0.020 


Charbonnel et al. (1999) 


Geneva 


0.691 


0.019 


D'Antona & Mazzitelli (1997) 


DM97 


0.699 


0.019 


Swenson et al. (1994) 


SW94 


0.703 


0.020 


Siess et al. (1997) 


SFD97 


0.703 


0.020 


this work 





at the base of the convective enveloppe : Afbasc 



0.981999 M , 
10 6 K and p ha 



Rba 



0.73322 R s , T b a 



1.999 x 



1.39610" 2 gcm- 



These numbers are in very close agreement with other 
recently published standard models (see e.g. Brun et al. 
1998, Morel et al. 1999, Bahcall & Pinsonneault 1996). 

In our grid computed with overshooting, we found that 
the 1 M Q model maintains a small convective core of ~ 
0.05M Q during the central H burning phase. This would 
consequently be the case for the fitted sun but present day 
hclioscismic observations cannot exclude this possibility 
(e.g. Provost et al. 2000). 

4. Comparison with other works 

In this section we compare our PMS tracks with the com- 
putations made available by the groups listed in Tabic |l|. 
These comparisons show the accuracy of the EOS and the 
pertinence of our models in the domain of VLMS. 

Pre-main sequence tracks differ from one group to an- 
other due to differences in the constitutive physics (EOS, 
convection), but also in the treatment of the surface 
boundary conditions. In the last decade, a tremendous 
amount of work has been done to better understand the 
physical processes taking place in the extreme regime of 
high density and low temperature encountered in the inte- 
rior of VLMS. Three main EOS have emerged: the MHD 
EOS used by the Geneva group and by DM97, the SCVH 
EOS (Saumon et al. 1995, hereafter SCVH95) used by 
Baraffe et al. (1998, hereafter BCAH) and the OPAL 
EOS (Rogers et al. 1996) used by DM97 in the regime 
of high temperature. Comparisons between these different 
EOS (SCVH95, Rogers et al. 1996, Chabrier & Baraffe 
1997, Charbonnel et al. 1999) showed a good agreement 
within their domain of validity. For example the MHD and 
SCVH EOS give very similar results in the mass range 
O.4£M£O.8M . 

Figure ^ shows that, the morphology of our tracks is 
very similar to the BCAH and Geneva models, for all the 
considered masses. This is a strong evidence that we cor- 
rectly follow the thermodynamics involved in these ob- 
jects. This is not surprising since, as mentioned in Sect. 
2.3.2, we treat the non-ideal effects through analytic fits to 
the MHD EOS. In comparison, our previous tracks (long 
dashed-dotted lines) were systematically too red and their 
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Fig. 3: HR diagrams comparing the different PMS tracks 
for the models listed in Table [l] 



paths were very different from the new ones. This is due in 
part to the fact that H2 was not accounted for in the com- 
putations and also because the treatment of the non-ideal 
effects was less accurate. As a consequence, the internal 
structure of our old models was less centrally condensed 
and the stellar radius slightly larger compared to our new 
models. 

At the lower mass end (M=0.10 and 0.20 M Q ), the sim- 
ilarity with BCAH models is striking, as effective temper- 
ature differences are 5s 50K (much less than observational 
uncertainties along the Hayashi tracks!). In particular, this 
shows that our boundary conditions (Sect. |2.2D are in rea- 
sonnably good agreement with the much more sophisti- 
cated atmosphere models used by BCAH. We see however 
that in the domain of very low mass stars (M ^ 0.4 M Q ) 
strong morphological differences are present. More specifi- 
cally, we observe (Fig. |) that the DM97 and SW94 tracks 
present a strong inflexion in their Hayashi path which is 
not reproduced by the other groups. Several reasons can 
be put forward to explain these differences. First DM97 
uses a convection model which differs substantially from 
the commonly used standard MLT prescription. As shown 
by D'antona and Mazzitelli (1998), this can result in T e g- 
differences as large as ~ 200K along the Hayashi line 
(see also below). Secondly, in this mass range the mor- 
phology of the evolutionary tracks is particularly sensi- 
tive to the EOS and to the different treatments of the 
thermodynamics. Finally, their use of a grey atmosphere 
approximation is certainly not satisfactory for modeling 
cool stars whose emergent spectrum diverges significantly 
from a black body. In order to estimate the impact of us- 
ing a grey atmosphere treatment in our stellar evolution 
code, we have also computed sonic additional tracks for 
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different masses. We report that the non-grey atmosphere 
models are systematically cooler than the grey ones, as 
a result of molecular blanketing in the outermost layers 
(e.g. Chabrier et al. 1996). This probably also explains 
why the Geneva models, which use the same treatment 
of convection with a similar value for MLT parameter 
(a = 1.6 ± 0.1), are systematically hotter than ours. 

For higher mass stars (M <; 0.6 M Q ), the morphol- 
ogy of the tracks becomes quite similar as the non-ideal 
physical effects and the molecular absorption in the at- 
mosphere vanish. The shifts in T e g between the different 
tracks mainly result from the different chemical composi- 
tions and/or values of a in the MLT and can be as large 
as 200K. In particular, the Hayashi tracks of the BCAH 
models become increasingly cooler as the stellar mass in- 
creases compared to all the other tracks. This most prob- 
ably results from their smaller value of the MLT param- 
eter (a = 1), the effect of which is only significant above 
~ 0.6M Q (Chabrier and Baraffe 1997) and is to shift the 
tracks towards lower T e g. The fact that species heavier 
than H and He are not accounted for in the BCAH EOS 
could be part of the explanation as well. Finally, for the 
most massive stars (M > 1.0 M ), the morphology of the 
tracks is very similar. We notice however that our models 
are slightly cooler than the DM97 tracks by ~ 100K. 

A comparison of the isochrones shows that our tracks 
are generally slightly more luminous than those of oth- 
ers. Consequently, for a given luminosity, the stellar age 
estimated from our models will be smaller, the relative 
effects being more pronounced for younger stars. The dis- 
crepancy between the different isochrones is particularly 
strong below 10 6 yr, where an age determination remains, 
in any case, very misleading. This is in part due to the 
fact that early in the evolution, the results still depends 
on the initial state of the star (Tout et al. 2000). 

Aside from the above-mentioned explanations for the 
discrepancies between the tracks, other physical ingredi- 
ents usually not reported in the literature may also af- 
fect the evolutionary path of a PMS stars. Consequently, 
the origin of some of the observed discrepancies cannot 
be clearly identified. To illustrate this purpose, we will 
now demonstrate that, even within the MLT, different pre- 
scriptions to solve the third polymonial equation for the 
convective gradient give rise to different degrees of supera- 
diabadicity and thus different location of the star in the 
HRD. Indeed, the T e g location of a fully convective PMS 
stars is very sensitive to the degree of superadiabadicity in 
the sub-photospheric layers which in turn depends on the 
prescription used to compute the convective temperature 
gradient (see e.g. D'antona & Mazzitelli 1998). We per- 
formed a series of tests using two different formalisms for 
the computation of the convective gradient, namely the 
Cox (1984, Chap. 14) and the Kippenhahn and Weigert 
(1991, Chap. 7) formalisms. The latter, used in our grids, 
leads in the superadiabatic region to a temperature gra- 
dient ~ 20% higher than the one coming out from the 
Cox formalism. This translates into effective temperature 



shifts of the Hayashi tracks of 100 to 300K, depending 
on the stellar mass. The Cox solutions have systemati- 
cally smaller radii and the Hayashi lines are consequently 
bluer. Finally with the receding of the convective enve- 
lope, the degree of superadiabadicity decreases and the 
temperature shift reduces to ~ 100K. These differences 
mainly result from the assumed "form factor" used in the 
MLT formalism. Indeed, the shape of the convective cells 
directly affects the radiative loss efficiency and a change 
in the form factor can be compensated by a modification 
to the a parameter. We refer the reader to Henyey et al. 
(1965) for a detailed discussion of the arbitrariness of some 
constants used in the MLT. Among other numerical tests, 
let us also mention that the evolution of the structure of 
PMS stars does not depend at all on the prescription used 
to write the gravothermal energy production rate e grav 
(i.e. using the internal energy, the entropy or the pres- 
sure; see Kippenhahn & Weigert 1991, Chap. 4). This is 
due to our completely consistent treatment of the thermo- 
dynamics with our EOS formalism. 

Given the relatively good agreement between the dif- 
ferent EOS involved in this comparison, it comes out that 
the physical and numerical treatment of convection is cer- 
tainly the most influent parameter for the modeling of 
fully convective PMS. It can account for T e g differences 
as large as 300K. Faced to our poor knowledge of con- 
vection, large discrepancies between the different sets of 
tracks are unfortunately inevitable. 

5. Internet interface to stellar model requests 



The internet site located at http:/ /www-laog.obs.ujf-greno- 



|ble.fr/activites/starevol/evol.html| offers several services to 
take advantage of our stellar model database. Among 
these facilities, it is possible to compute an isochrone, to 
determine the stellar parameters of a star knowing its po- 
sition in the HR diagram, to draw a specific HR diagram 
and of course, to retrieve all the evolutionary files. 

For the computation of isochrones, the user specifies 
the metallicity, ages and selects the mass tracks enter- 
ing the computation of the isochrone. The results are dis- 
played on the screen and can be saved in a file. We also 
give the luminosity, effective temperature, radius as well as 
the colors and magnitudes corresponding to the location 
of a given star at any specified age. We use the conversion 
table provided by Kenyon and Hartmann (1995) and dis- 
play the different colors in the Cousin system. It is also 
possible to ask for further information and get, for exam- 
ple, the surface Li abundance, the central temperature or 
the extent of the convective envelope for the set of stars 
selected to compute the isochrones. 

The other facility is the determination of stellar pa- 
rameters. The user specifies the coordinates of a star in 
the HR diagram using either the luminosity or magnitude 
and T e g of colors and then, the program computes the the- 
oretical track that passes through this observational point. 
This procedure allows the determination of the stellar age, 
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mass and radius assuming the star is on its pre-main se- 
quence track. If requested, additional information such as 
surface chemical composition, central properties or mo- 
ments of inertia can also be provided. 

The third facility deals with graphic outputs. In this 
page, the user can plot HR diagrams in any combination 
of color and magnitude. It is also possible to zoom in by 
stating limits, to superpose isochrones and/or the ZAMS. 

6. Conclusion 

We presented new pre-main sequence evolutionary tracks 
for low- and intermediate-mass stars. Comparisons of our 
models with other PMS tracks indicate rather strong dis- 
crepancies in the regime of very low-mass stars {M ^ 
0.5 Mq). Our tracks are very similar, in morphology and 
effective temperature, to those computed by Baraffe et al. 
(1998) and the Geneva group (Charbonnel et al, 1999). 
This similarity is a strong indicator that our treatment 
of the EOS and boundary conditions are correct. Con- 
versely, we report strong deviations of these tracks with 
the D'Antona and Mazzitelli (1997) ones, especially be- 
low 0.3 Mq. For higher mass stars, the morphology of 
the tracks are similar and effective temperature differ- 
ences Ss 200K are noted, partly due to different mixing 
length parameter values and chemical compositions. Com- 
parisons of the isochrones indicate age determination re- 
mains, in any case, very uncertain below 10 6 yr. Finally, 
we present our internet server which provides several fa- 
cilities to use and take advantage of our large database of 
PMS stellar models. 
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